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We discuss the role of the proximity force approximation in deriving limits to the existence of 
Yukawian forces - predicted in the submillimeter range by many unification models - from Casimir 
force experiments using the sphere-plane geometry. Two forms of this approximation are discussed, 
the first used in most analysis of the residuals from the Casimir force experiments performed so 
far, the second recently discussed in this context in R. Decca et at, [Phys. Rev. D 79, 124021 
(2009)]. We show that that the former form of the proximity force approximation overestimates the 
expected Yukawa force and that the relative deviation from the exact Yukawa force is of the same 
order of magnitude, in the realistic experimental settings, as the relative deviation expected between 
the exact Casimir force and the Casimir force evaluated in the proximity force approximation. 
This implies both a systematic shift making the actual limits to the Yukawa force weaker than 
claimed so far, and a degree of uncertainty in the a — X plane related to the handling of the various 
approximations used in the theory for both the Casimir and the Yukawa forces. We further argue 
that the recently discussed form for the proximity force approximation is equivalent, for a geometry 
made of a generic object interacting with an infinite planar slab, to the usual exact integration of 
any additive two-body interaction, without any need to invoke approximation schemes. If the planar 
slab is of finite size, an additional source of systematic error arises due to the breaking of the planar 
translational invariance of the system, and we finally discuss to what extent this may affect limits 
obtained on power-law and Yukawa forces. 

PACS numbers: 04.50.-h, 04.80.Cc, 03.70.-|-k, 12.20.Fv 



I. INTRODUCTION 



Several unification schemes merging gravity and the standard model of strong and electroweak interactions predict 
the existence of short-range forces with coupling strength of the order of Newtonian gravity [ij . Efforts to evidence 
a Gfth force have been envisaged regardless of any concrete unification scheme since various decades 0, Q , and there 
are compelling reasons to improve our limits especially in the largely unexplored submillimeter range. Constraints 
in both coupling and range for these interactions have been obtained with various experimental setups, including 
the recent configurations using a disk-shaped torsional balance parallel to a rotating flat surface [1, H, @, 0, or 
micromcchanical resonators in a parallel plane geometry [^, [lO, [HI, [13, [ll] • Due to the surge of activity in the 
study of Casimir forces, limits have been also given in the submicrometer range based on the level of accuracy between 
Casimir theory and experiment. However, unlike the case of experiments performed between bodies kept at larger 
distances, the use of the parallel plane geometry on such small lengthscale has been proven to be challenging in terms 
of parallelism [ij, [Hi [I3| , and therefore the attention has been focused on the analysis of the residuals in the Casimir 
theory-experiment comparison involving the sphere-plane configuration. 

Dedicated efforts to obtain limits from sphere-plane Casimir experiments have involved the use of the so-called 
Proximity Force Approximation (PFA) [VJi, [l^, [13, [Ml, H^, H^, which allows to map the force Fsp between a 
sphere of radius R and a plane located at a distance a from the sphere into the energy per unit area Epp of the 
parallel plate configuration, namely Fsp{a) = 2nREpp{a) [1^. This approximation is believed to be valid in the limit 
a <C i? and to hold with a high degree of accurac y fo r forces between entities concentrated on the surfaces, such as 
electrostatic or Casimir forces between conductors |25l . [26l . [27| . Obviously, in order to test how well PFA approximates 
the exact force, one needs either to compute the interaction exactly or at least to assess reliable bounds. For the 
electrostatic sphere-plane interaction, the exact analytical result for the force is well-known and has a closed form [2^ , 
such that deviations from PFA can be readily anal yze d. For the Casimir sphere-plane interaction, the exact force has 
been computed only very recently, both for ideal pol . Isol . Isij and real metallic plates [s^l ■ Available analytical and 
numerical results seem to indicate that, at least for zero temperature and within the used plasma model, deviations 
from PFA applied to the sphere-plane Casimir interaction are small, of the order of 0.1% or higher, in recent Casimir 
experiments aiming to put limits to Yukawa interactions. 
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It has been argued in [33| that the appHcation of the PFA to forces acting between entities embedded in volumetric 
manifolds, such as gravitational forces or their putative short-range components, is in general invalid and has to be 
carefully scrutinized in each specific configuration. Based on this suggestion, a recent reanalysis of the PFA in the 
case of gravitational and Yukawian forces has been discussed in [33 | . The main conclusion of this reanalysis is that "a 
confusion with different formulations of the PFA" existed in the previous literature, and that "care is required in the 
application of the PFA to gravitational forces" . This confusion is stated to originate from a specific form of the PFA 
used so far, to be contrasted with a more general formulation of the PFA. In [34| it is also claimed that the difference 
between the two PFAs is negligible in the actual configuration used to give the allegedly best limits obtained in the 
100 nm range 

In this paper we further discuss the meaning of the PFA in the case of volumetric forces. We argue that the 
discrepancy between the two forms of the PFA is a significant source of error in the determination of bounds on 
parameters of Yukawian forces from force residuals in Casimir sphere-plane experiments performed so far that used 
PFA to model such non-Newtonian forces. We then show that the general form for the PFA discussed in [s^l is simply 
a different choice of the infinitesimal volume for integrating the force due to an extended object, and coincides with 
the exact result only in the case when one of the two surfaces is an infinite planar slab (or semispace). The level of 
approximation in using the two PFAs for Yukawa forces in the sphere-plane geometry is of the same order of magnitude 
as the Casimir theory-experiment comparison, that uses PFA to compute the sphere-plane Casimir force (as already 
noticed in 3J|). Therefore, since such a comparison provides force residuals that are in turn compared against the 
theory of Yukawa forces to obtain limits on its a — A parameter space, the use of these subsequent PFA approximations 
of comparable level of approximation provides a possible source of systematic error, not carefully accounted for so far. 
We also argue that other volumetric effects not directly related to the PFA, such as the finite size of the planar surface 
used in the actual experiments, may provide a source of systematic error not taken into account so far, which strongly 
affects the limits to power-law forces [H,|35|, but should not be a major source of concern on limits to Yukawa forces. 
We believe that, considering the various number of complications related to the sphere^lane geometry, upgraded 
versions of parallel plates experiments such as the ones discussed in d, [^, [l3, [HI, [13, [iS] could provide limits on 
Yukawian and power-law forces in the submicrometer range more immune to a set of systematic errors characteristic 
of the sphere-plane configuration. 



II. PROXIMITY FORCE APPROXIMATIONS AND VOLUMETRIC FORCES BETWEEN EXTENDED 

OBJECTS 

In order to introduce the notation and as a prelude to our discussion, we briefly summarize the results contained 
in (s^ l . The actual experimental configuration used in [2]| is not a parallel plate geometry, rather it is a sphere-plane 
geometry, and the PFA is used to map the force between a sphere and a plane Fgp into the energy per unit area of 
the parallel plate configuration Epp 

F,p{a)^2TTREpp{a), (1) 



where R = y/R^Ry is the geometrical average of the principal radii of curvature of the spherical surface evaluated at 
its point of minimum distance from the plane. In the experiment reported in plj . the force is measured by looking 
at the frequency shift of a mechanical resonator, as customary in atomic force microscopy [1^, and as first reported 
in the context of Casimir force measurements in (37j . The frequency shift is proportional to the gradient of the force, 
and therefore 

Azy2 = ^ ^-^^P ^ ^ ^^PP ^ ^ p (2) 
An'^m da 27rm da 2i:ni 

where Ppp is the plane-plane pressure, and m is the mass of the resonator. The measure of the frequency shift can then 
be mapped, via use of Eq. ([T]), into the equivalent pressure exerted between two fictitious parallel plates mimicking 
the actual sphere-plane geometry. Within the validity of Eq. ([T|), this is a valid assumption for the case of forces 
acting between surfaces, such as electrostatic forces between conductors or Casimir forces. 

A first sign of the fact that there can be issues with the PFA in dealing with volumetric forces, such as the 
hypothetical Yukawian forces of gravitational origin, is manifested by noticing that the exact formula for the Yukawa 
force between two infinite parallel slabs depends on the thicknesses of both slabs, which implies that the PFA formula 
applied to the volumetric Yukawa force in the sphere-slab configuration also depends on both thicknesses (and on the 
sphere radius). However, the exact sphere-slab force obviously depends only on the slab thickness and on the sphere 
radius - it does not, and cannot, depend on the thickness of the metaphysical slab introduced in the virtual mapping 
to the parallel geometry. Indeed, consider the Yukawa potential energy for two pointlike masses mi and m2, located 
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FIG. 1; Schematics of the integration for the sphere-slab configuration according to the usual slicing along horizontal infinites- 
imal slabs used for the exact calculation (left), and slicing using vertical columns (right) as used in the EPFA calculation. As 
far as the surface facing the sphere is an infinite plane (and therefore translational invariance of the potential V due to the 
plane is satisfied), these merely correspond to two different and equivalent choices for the integration volume. 

at positions ri and r2 respectively, 

g-|r2-ri[/A 

C^Yu(ri,r2) = ~aGmim2 -. p, (3) 

|r2 -ril 

where, as usual, the strength of the Yukawa interaction is parameterized in terms of Newton's gravitational constant 
G through a dimensionless quantity a, and A is its range. Assuming that the Yukawa interaction is additive, once 
integrated over two infinite, homogeneous parallel slabs separated by a distance a, one derives the corresponding 
pressure Pyu, 

PYu(a) = -27TaGpip2X^e--/\l - e-^i/^)(l - e"^-/^), (4) 

where Di and D2 indicate the thickness of each slab, pi and p2 their densities. The exact Yukawa interaction in the 
sphere-slab geometry can be readily computed, assuming additivity. The result is [s^ 

F^^-='(a) = -47r2aGpip2A3i?e-"/^(l - e-°'^^){l - \/R + e~'^^'^ + e~^'^'^\/R). (5) 

As we mentioned above, most recent experimental works on limits to extra-gravitational forces from sphere-plane 
Casimir measurements used the usual PFA. In this approximation, the Yukawa force between a homogeneous sphere 
and an infinite homogeneous slab of thickness Di is 

F^^^ia) = 27ri?PYu(a) = -ATi''aGpiP2\^ Re-^-'^il - e'^^'^){l - e,-^^/^). (6) 

In this case one needed to consider, in order to map the actual sphere-plane configuration into a parallel plate 
geometry, a fictitious upper plate of thickness D2 large enough, i.e. much larger than the explored Yukawa range 
{D2 S> A). Again, this situation may appear disturbing to whoever believes that any experiment-theory comparison 
should not rely on the introduction of arbitrary parameters not having a tangible, measurable counterpart in the 
concrete experimental setup. Clearly the PFA prediction Eq. ^ fails to give the exact result Eq. ([5]) in the range of 
its supposed validity, a R, and it is necessary to assume X R, D2 in order for PFA to tend to the exact result. 
But in this limit the volumetric nature of the interaction is lost, since the atoms in the "bulk" no longer contribute 
appreciably to the total force. Likewise, PFA fails to give the exact Newtonian interaction between the sphere and 
the slab, even in the range of its supposed validity a <^ R. 

The authors of [s^l consider the most general formulation of PFA [11], that we will call "exact" PFA formula 
(EPFA) to distinguish it from the usual PFA approximation, described above. In the EPFA the force between two 
compact bodies is expressed as the sum of forces between plane parallel surface elements dxdy. The z component of 
the force is 



i^™(a) = j jjxdy P{x,y,z{x,y)), 



(7) 
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where P{x,y, z{x,y)) is the pressure between two parallel plates at a local distance z{x,y) = Z2{x,y) — zi{x,y) > 
{zi{x,y) are the surfaces of the two bodies), a is the distance between thcni (smallest value of z{x,y)), and a is the 
part of the (x, ?/)-plane where both surfaces are defined (see Fig. 1). The EPFA prediction for the Yukawa force 
between a sphere and an infinite planar slab is [33 | 

i^™(a) = -47r^aGpip2X''Re-/\l - e-^i/^)(l - X/R + e'^^/^ + e'^^AA/i?), (8) 
which coincides with the exact result of Eq.(l5]). In Fig. 2 we plot the ratio r] 

_ f|PF^ _ l-X/R + e-^'^/^ + e-^^'/^X/R 

as a function of the Yukawa parameter A, for different values of the sphere radius keeping fixed D2 oo (left plot), 
and for different values of D2 keeping fixed the sphere radius at i? = 150/im (right plot). Note that rj is independent 
of the sphere-slab separation a. When D2 ^ A, as surely realized in the left plot, PFA always overestimates the EPFA 
result, i.e. 77 < 1 (similarly to how the PFA overestimates the exact Casimir force in the sphere-plane geometry). Note 
also that when D2 ^ X the atoms in the "bulk" of the two bodies do not contribute appreciably to the Yukawa force, 
thereby making it effectively of a surface character (i.e., non- volumetric) , as in the case of Casimir or electrostatic 
forces. Instead, for values of D2 — lOA (or smaller) the volumetric character of the Yukawa interaction is manifest, 
and T] is no longer less than one (right plot) . In this case the PFA applied to Yukawa forces is invalid and in the data 
analysis one should at least declare the value of D2 imagined for which the limits are assessed. Overestimating the 
Yukawa force leads to stronger limits for the coupling constant a for a given A with respect to the proper use of the 
EPFA. Considering the relatively small margins of improvement reported recently (see for instance Fig. 3 in [2l[) a 
systematic shift due to the use of the PFA instead of the EPFA may lead to significant changes for the exclusion region 
in the a — A plane. Both plots show that the use of PFA instead of EPFA is unreliable especially in the region near 
or above A=100 nm. Unfortunately the region in between 100 nm and few fim is also the one directly explored with 
the Casimir force experiments, since actual measurements take place in this range of distance between the involved 
objects. It is known that the best limits on Yukawa interactions can be set for A of the order of the actual explored 
distance between the two bodies ~ a, and the extrapolation of the measurements to smaller A is affected by the fast 
growth of the bounds as a oc exp(a/A). 

The fact that the EPFA Eq. ^ gives the correct exact results for the Yukawa (and also gravitational) force in the 
sphere- plane configuration is in fact a trivial consequence of the additivity of these interactions and of the translational 
invariance of the infinite plane surface, the shape of second surface (in this case a sphere) being irrelevant. Indeed, 
the EPFA is just a different parameterization of the exact formula of addition of forces between particles. On one 
hand, the exact interaction energy between a test mass and an infinite slab (or half-space) depends only on the normal 
coordinate z, being independent of the in-plane coordinates x, y by symmetry. Hence, the potential due to the infinite 
slab is V(x,y,z) = V{z). For a body {e.g. a sphere) of mass density p{x,y,z), additivity implies that the total 
interaction energy can be obtained as 

C/body = J dxdydz p{x, y, z) V[z), (10) 

and similarly for the force. Since V depends only on z, it is convenient to compute the integral by adding forces at 
different slices at constant z, i.e. considering infinitesimal slices in z, then evaluating first the potential energy of a 
slice of the body parallel to the plane at a distance z 



W{z)^ j dxdy p{x,y,z)Viz), (11) 

and then integrating along z the quantity W(z) one obtains the exact expression for the body-plane interaction ?7body 
(see Fig. 1 left). On the other hand, the EPFA states that the interaction energy between the body and the plane 
is obtained from slicing the body into cylinders perpendicular to the plane (see Fig. 1 right) and integrating the 
cylinder-plane interaction energy along the portion of the body that faces the plane (i.e., one must integrate over the 
surface a on the plane that is the normal shadow of the body). The potential energy of this column of the body 
centered around {x,y) is 

G{x,y) = J dzp{x,y,z) Viz), (12) 

and then integrating along x, y the quantity G{x, y) one gets the EPFA expression for the body-plane interaction 
t^body (sec Fig. 1 right). Again, since the interaction is additive and it does not depend on the x,y coordinates, this 




FIG. 2: (Color online) Comparison between the EPFA and the PFA versus the range of the Yukawian force for a homogeneous 
sphere above an infinite homogeneous slab. (Left) Plot of 77 versus A in the case of spheres of radius R =150, 100, and 50/im, 
in the limit of D2 — > 00. (Right) Same plot but for a sphere of radius 150/im and finite values of the metaphysical parameter 
D2 representing the thickness of the upper slab introduced in [sj for the comparison of the actual sphere-plane geometry to 
the parallel plate case using the PFA. 



integral is exactly equal to the previous one: we are simply integrating the same function U{z) over the body using a 
different parameterization of the volumetric integral. The same holds for the force between any body (not necessarily 
a sphere) and an infinite slab (see also [3§|). 

However, when none of the two bodies is an infinite slab (e.g., two spheres of radii Ri, mass densities pi{x,y,z), 
separated by a distance a along the z direction), translational invariance along a; — y is obviously broken, and the 
EPFA does not coincide with the exact formula. The exact interaction energy U{x,y,z) between a source body and 
a test mass at position {x, y, z) can be easily computed. For instance, for the spherical source body U depends only 
on r = (a;^ + + z^Y^^, with the origin of coordinates at the center of the sphere. Integrating U{x,y, z) over the 
volume of the second body one gets the exact result. For example, for the two-spheres case the gravitational energy 
depends only on the center-to-center distance and scales as 1/a. Let us compare this known exact result with the 
EPFA prediction. One slices each body in cylindrical slabs, calculates the slab-slab interaction energy Uss{z) that 
depends on the local distance z between the slabs, and finally one adds up these contributions over the shadow of one 
of the bodies on the other one. It is clear that EPFA cannot give the exact result since Uss is translational invariant 
but the exact U is not, and EPFA fails to predict the exact 1/a dependency. Therefore the EPFA formula for the 
energy and force of additive two-body interactions is a trivial reparanietcrization of the exact result when one of the 
bodies is an infinite plane (or slab). For other geometries, EPFA fails to give the correct result as a consequence of the 
broken translational invariance. In particular, this is the case of a sphere above a finite size slab, like in experiments, 
especially those involving slabs of typical sizes comparable to those of the sphere. 

In the experimental configuration used in [2lj various substrates are present on the sphere and on the slab. Imagining 
that the layered slab is infinitely long, the Yukawa potential at a distance z from the top layer due to the slab is 

y^Jz) = -27r«GA2e-^A L/'g-A'//A(gA'//A _ 1) + g-(A'/+A;)/A(gAi/A _ 

^^g-(A'/+A; + Di)/A(gDi/A _ _ (^3) 

Here A" and are the thickness and density of the top layer, A[ and p[ are the thickness and density of the middle 
layer, and Di and pi and the thickness and density of the lower part of the layered slab. The last factor in Eq. p^ 
can be considered a sort of effective density of the planar surface, in which the various densities are weighted by their 
thicknesses in units of A (indeed yielding their arithmetic average in the case of A', A", Di <^ X). 

We can compute the exact expression for the Yukawa interaction energy between the layered infinite slab and a 
layered sphere of mass density p2{x,y,z) using Eq. ([TO)) . As discussed above, this exact computation will trivially 
coincide with the EPFA expression. Let R and p2 be the radius and density of the sphere, A'2 and p'2 the width and 
density of the inner layer on the sphere, Aj and the width and density of the outer layer, and a the distance from 
the outer layer of the sphere to the top of the layered slab. The total EPFA Yukawa interaction energy can be written 
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as a sum of contributions from each layer on the sphere, Uy^ 

rR 



A.EPFA 



U^^, where 
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r^dr / desmOp'^V^'^iz) 



Note that, instead of using horizontal or vertical slicings for the volume integration as done in Fig. 1 for the non- 



layered case, we use spherical slicings more appropriate for the layered sphere case. Here z 



Af. + i? - r cos 



denotes the vertical position of any infinitesimal mass element inside the layered sphere. Computing these integrals 
we obtain the EPFA expression for the Yukawa interaction energy between the layered infinite slab and the layered 
sphere 
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(14) 
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Note that when 



The EPFA expression for the corresponding force is Fy^ — —ui^y^ 
there are no layers on the slab {A[ = A" = 0) and no layers on the sphere (A2 = Aj = 0), then the expression for the 
force that follows from Eq. is identical to Eq. ([S]). On the other hand, the PFA expression for the force between 
the layered infinite slab and the layered sphere is Fy^ = 2TrRPy^, where Py^ is the pressure between two parallel 
layered slab, one identical to the previous slab, and a metaphysical slab of width D2 and density p2, covered by two 
layers of widths and densities identical to the ones of the layered sphere above. Using Eq. (|4]) for the various pairs of 
layers in the different slabs, we calculate the PFA expression for the layered sphere-slab force 



i;,A,PFA, 

^Yu (« 



(a) = -An^aGX^Re-"^^ x 



Pi 



=.-A'i'A 



P2 

"^^g-(A'/+Ai+Ai')/A(i _ e-^2/A) ^ ^^g-(A'/+Ai')A(i _ e^^'^^) + p^'e-^"/^(l - g-^^'/A) 



(15) 



To assess the effect of the multilayered structures, we have evaluated the ratio ija/v with ija = Py^^^^^ / Pyxf^^ s-J^id 
77 = Fy^^-^ / Fy^-^ , as a function of A for three radii of curvatures of the sphere, assuming for the PFA calculation a 
value of the metaphysical parameter D2 = lO^^m (sec Fig. 3 left). The effect of multilayers is to slightly flatten ?/a 
as compared to ?/ in the homogeneous case (Fig. 2 left). The dependence of the same ratio for a fixed value of R 
and different values of the metaphysical parameter D2 is shown in Fig. 3 right. Note that 77A is independent of the 
sphere-slab separation, just as 77 is. 

As discussed in the Introduction, the PFA used in all recent sphere-plane Casimir experiments for the Casimir 
theory-experiment comparison is expected to approximate the exact Casimir force within 0.1 %. This expectation 
comes from recent analytical approaches to the sphere-plane Casimir interaction [2^, [13, HH, [13] that, although formally 
exact, require the evaluation of the determinant of an infinite-dimensional matrix, which becomes a numerically 
demanding task, especially in the PFA regime, a <^ R, where larger and larger matrices are needed for convergence. 
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FIG. 3: (Color online) Ratio r\i\lr\ for the comparison of the EPFA and PFA schemes for the multilayered and corresponding 
homogeneous situation (obtainable by using Eqs. (14) and (15) with A2 = A2 = 0, P2 = p'i ~ 0, and replacing R with 
R + A'2 + A2 )• (Left) Ratio t^a/t? versus the range of the Yukawian force for different values of the radius of the inner sphere. 
The parameters for the layered sphere are A2 = 10 nm, A2 = 180 nm, p2=4.1 g/cm"', p2=7.14 g/cm^, and p2=19.28 g/cm"^. 



=7.14 



and 



The parameters for the layered slab are Di = 3.5/im, A'l = 10 nm, A" = 210 nm, pi=2.33 g/cm'^, p'i = 
pi=19.28 g/cm"^. In the evaluation of the PFA force Fy^^^^, a value of the metaphysical parameter D2 — 10*/im is used. 
(Right) Ratio 77A/7? versus the range of the Yukawian force for different values of the metaphysical parameter D2 and a radius 
of curvature oi R = IbOfim. 



Numerical computations of the exact, zero-temperature sphere-plane Casimir force using parameters for metallic 
spheres {R = lO/im and optical response modeled by the simple plasma model with plasma wavelength Ap = 136nm) 
show that deviations from PFA can be as large as 20% for the smallest a/i? w 0.5 studied numerically (see Fig. 2 of 
[32I). An extrapolation to smaller values of a/R using a cubic polynomial fit of the numerical data is also provided 
in j32l |. Assuming one can use it for the recent Casimir sphere-plane experiment [2l] | (with a radius of curvature 
R = 151. 3^m), gives a deviation from PFA of the order of 0.1% at the smallest value of a/R « 0.001 reached in 
the experiment (amin ~ 160nm). Since the limits to non-Newtonian forces are obtained using the residuals in the 
Casimir theory-experiment comparison, in order to meaningfully replace the exact formula of the Yukawa force with 
its PFA approximation, the level of accuracy between these two should be therefore a small fraction, for instance 10 
%, of the accuracy with which the Casimir force is controlled by using PFA rather than the exact expression for the 
sphere-plane Casimir force. If this condition is not fulfilled, the derived limits could be off also by a large, order of 
100%, correction. However, targeting a 10% accuracy level with respect to the Casimir theory-experiment accuracy 
implies deviations from 77 = 1 of 0.01%, which can be obtained, as seen in Fig. 2, only in the range of A below 100 
nm. The presence of substrates with different densities tends to mitigate the discrepancy between the EPFA and the 
PFA, as seen by the curves in Fig. 3, but there is an irreducible systematic factor even at small A. Indeed, in the 
limit A — ^ 0, we have 77A ~ 1 + (Aj + Ag)/^?, that, in the case of the experiment reported in [2l|, is equal to 1.00126, 
i.e. a correction already equal to 0.126 %. 

All these systematic sources of uncertainty could be even larger in experiments for which the radius of curvature of 
the sphere is not adequately optimized. Indeed, the use of spheres with smaller radius of curvature is affected more 
by this effect, as emphasized in the left plot of Fig. 2 and in Fig. 3 for the cases of i? = 50^m and lOO^m. Moreover, 
for small spheres the PFA approximation to the Casimir force itself is less accurate. The use of spheres with large 
radius of curvature is beneficial to reduce these sources of error in the experiment-theory comparison, but may face 
experimental issues recently identified in [40| and interpreted as due to deviations from an ideal spherical geometry 
(as proposed in [4l|) and/or a consequence of larger sensitivity to electrostatic patch effects (33ll4g|. 



III. BREAKING THE X-Y TRANSLATIONAL INVARIANCE 



In the previous section we have seen that translational invariance is crucial to make the EPFA reproduce the exact 
result, but in actual experiments such an invariance is obviously satisfied only approximately, leading to an additional 
source of systematic error related to the finite size of the surfaces, as we discuss here for both power-law forces, and 
for Yukawa forces. Instead of computing the more involved problem of the gravitational force between a sphere and a 
finite-size slab, wc consider here the simpler case of the gravitational force acting on a point-like test mass m2 above 
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FIG. 4: (Color online) (Left) Plot of the ratio ^at between the force F(r) oc exerted at the closest point of the sphere from 
the disk and the force exerted at the farthest point, versus the radius of the disk Rd constituting the planar surface of finite 
size, for four different exponents A*'. We assume a radius of the sphere equal to _R = 150/im, a sphere-plane distance of a=100 
nm, and that the center of the sphere is right above the center of the disk. Only the case of A'' = 2 (Newtonian gravitation) 
gives a ratio of unity in the large Rd limit. (Right) Plot of the ratio ^jv versus the exponent of the power force law A'' for various 
radii of the disk Rd- The cross indicates the ratio ^ = 1 for the case of an exponent N — 2 and an infinite plane; it is provided 
as a help to the eye to better show the convergence in the case of A'' = 2 of ^jv to unity with disks of progressively larger radii. 



the center of a disk of thickness Di, radius Rd, and mass density pi. We obtain 

rO rRd 



Fg{z2) 



-27rGpi?T!,2 



dzi 
Di Jo 



rdr- 



Z2 - Zl 



+ (Z2 - zi)2]3/2 



(16) 



where Z2 is the distance between the test mass and the disk. The force becomes independent of Rd only in the Umit 
Rd 3> Di, Z2 (in which case it is also independent of Z2)- In order to assess the different forces acting on the various 
parts of a sphere in the presence of a disk of finite radius we evaluate the ratio between the forces exerted at the point 
of the sphere closest to the plane {z2 = a) and the farthest point (z2 = a + 2R). This quantity is simple to evaluate 
yet provides a practical figure of merit for how much the extended geometry of the sphere is affected by the finite size 
of the disk. This gives a ratio S^g = Fg{z2 = a)/ Fg{z2 = a + 2R): 



211/2 



f3 + [(2 + 7)2 + k2]1/2 _ [(2 + 7 + /?)2 + k2]1/2 



(17) 



where we have defined P = Di/ R, •y = a/R, and k = Rd/R- 

This is a large correction, of the order of 300%, if a disk of radius equal to twice the radius of the sphere {Rd = 2i?), 
in a geometrical setting not dissimilar from the one used in [2ll |. is considered. Since the experiment in f2l| is anyway 
insensitive to the gravitational force, like any experiment performed in the micrometer range, this is not a major 
practical concern. However, in the case of a more generic power law such as = —Kpim2/r^ we get: 
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apart from the cases of A'^ = 1 and = 3 in which logarithmic integrations occur. In these two cases one obtains 
■KKpim2 



(18) 
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which in the limit Rd 



oo becomes independent of z, and 
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FIG. 5: (Color online) Plot of the ratio ^yu between the Yukawa force exerted at the closest point of the sphere from the disk 
and the force exerted at the farthest point, versus the radius of the disk Rd constituting the planar surface of finite size, for 
three different values of the Yukawa range A. As before, we assume a radius of the sphere equal to i? = 150^m, a sphere-plane 
distance of a=100 nm, and the center of the sphere right above the center of the disk. The values of A are chosen to be 100 ^m, 
500 fim, and 1000 /im. In the last case the force may be considered as a long range one and the farthest point on the sphere is 
also contributing almost as the closest one. For A progressively smaller than the radius of the sphere the ratio ^yu gets larger 
and larger and the dependence on Rd is not appreciable. 



which in the hmit Rd ^ oo behaves as hi(l + D\/ z^). Notice that the fact that the force is independent on the 
distance from an infinite plane is only characteristic of forces scaling with the inverse square of the distance, such 
as the gravitational force, making the integration of the force trivially geometrical. In the general case n 2, even 
in the situation of a sphere in front of an infinite plane, different points of the sphere will feel different forces, with 
the farthest point feeling smaller (larger) force for a power law exponent larger (smaller) than 2, as a consequence 
of the interplay between the solid angle and the distance scaling of the force, which makes peculiar the N = 2 case 
as expressed by the Gauss law. This is shown in Fig. 4 (left) for the cases of = 1,2,3, and 4, with the ratio 
£,N between the forces evaluated at the top and at the bottom of the sphere, and in Fig. 4 (right) by showing the 
same ratio versus the power law exponent for different values of the radius of the disk. Therefore, when considering 
power-law forces as the ones discussed for instance around Eq. 2 in [2^, one should then take carefully into account 
the finite size of the plane in deriving limits to these forces [431 . 

Finally, we discuss the effect of the finite size of the planar surface in the case of Yukawa forces. The potential 
energy of a pointlike particle of mass m2 located at height z along the axis of a planar disk surface of radius R^, 
density pi, and thickness Di, is 



Uy^iz) 



-aGpim2 /"^ dr r f [(z2 - z^f + ^2]i/2g-V(-2-.i)^+rVA 

Jo J-Di 



(21) 



and the related force is 
Fyvl[z) = - 



-27raGpiTO2Ae-^'/^(l - e-°^l^) 



+ 



i?rf[l + (z/i?rf)2]i/2(l^e-^iA) Rd [l + ((z + i?i)/i?d)2]i/2(i_e-^iA) 



(22) 



where the finite size terms appear as corrections to the indefinite plane formula originating by the first term alone. 
As before, we introduce as figure of merit the ratio ^yu = ^Yu(-22 = a)/^Yu(22 = a + 2i?). This ratio is very large 
for realistic configurations, expressing the short-range nature of the force. Indeed, even in the infinite plane limit we 
have a ratio of ^yu = e^^^^ ~ gSOOO ^j^g ^^^^ sphere of radius R = 150/im at a A = 0.1//m. The dependence 
on the disk radius become significant only at values of A comparable to the radius of the sphere, as shown in Fig. 5. 
The presence of suppression factors for the farthest point of the form e^^/^ makes very insensitive the Yukawa force 
to the finite size of the disk (for previous considerations, sec also (isj). 
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IV. CONCLUSIONS 

Our analysis, although confirming some of the results already discussed in [s^ . draws quite different conclusions 
from the common outcome. In particular, we argue that the application of PFA to volumetric forces is not rigorous 
if considered in its original formulation applied so far to compute th e sp here-plane Yukawa interactions in the most 
sphere- plane Casimir force experiments reported in [l^, [I^, [2O, [HI, [13, [111- ap plic ation to volumetric forces is 
instead of trivial nature if considered in the exact formulation EFPA discussed in [3J], since the latter is identical 
to the exact calculation, just differing in the choice of the infinitesimal integration volume. We have shown that the 
usual PFA is an invalid approximation to compute volumetric forces. In particular, it does not reproduce in its usual 
range of validity exact known expressions for gravitational and Yukawa interactions in non-translational invariant 
geometries, such as sphere-finite size slab or sphere-sphere configurations. The "exact" PFA is the exact expression 
for any additive two-body interaction when one of the bodies is translational invariant, as is the case for a sphere 
in front of an infinite homogeneous slab or half-space. For non-translational invariant geometries, EPFA also fails to 
give the exact result for volumetric interactions, even in the regime of parameters where it is assumed to be valid, 
therefore also being an invalid approximation for volumetric forces. 

The difference between the two formulations of the PFA is shown to affect significantly the limits obtained so far 
unless one considers a regime of Yukawa range so small, A <C i?, £'2, that the approximation of a surface force {i.e. 
neglecting the Yukawa force due to the atoms in the "bulk" of the two bodies, therefore manifestly of non- volumetric 
character) holds. By using the PFA the Yukawa force is overestimated and therefore the limits in the a — A plane 
become more stringent than by using the exact force estimated via EPFA. Moreover, the use of PFA instead of 
the EPFA for the parameters of the experiment supposed to provide the strongest limits to Yukawian interactions 
[2l| occurs with an accuracy of the same order of magnitude with which the exact Casimir force is expected to be 
also approximated by the corresponding PFA. On one hand, since the Casimir theory-experiment comparison provides 
force residuals that are in turn compared against the theory of Yukawa forces to obtain limits on them, the use of these 
subsequent PFA approximations of comparable level of approximation provides a possible source of systematic error, 
not carefully accounted for so far. On the other hand, both Casimir and Yukawa PFAs overestimate the respective 
exact forces, and therefore the systematic source of error introduced by their use might be less critical than expected 
on first principles. In any case, it is important to assess as much as possible both sources of errors or, in alternative, 
to use exact expressions for the Yukawa and Casimir forces. Furthermore, a systematic source of error present even 
in the EPFA scheme for finite-size planar surfaces has been discussed and shown to be significant only for power-law 
forces. Our analysis suggests that future limits (or reanalysis of experiments already performed) on Yukawian forces 
should rely upon the use of the exact expression for the Yukawa force, as performed in (isl. |4^. 

It is also worth to point out that all the above considerations hold provided that the simple scenario of additive 
forces is assumed, which is valid in general for weak forces among atoms in the low density limit, such that correlations 
leading to fluctuating forces are negligible. However, the hypothetical Yukawa forces should be located in a regime of 
coupling constants intermediate between the gravitational (additive) force and the Casimir (non additive) force. It 
is not understood a priori if the Yukawa force is weak enough to make the additivity assumption reliable, and this 
should be kept in mind in future broad-range searches of these forces. 

Finally, considering the complications emerging in the sphere-plane geometry due to the presence of previously 
unidentified systematics such as the sensitivity to deviations from the ideal spherical geometry [ssl . |40| . [ill . l45j and 
possible effects of variability of the contact potential with distance [i^, [1^, [43, Tisl . l49j , it may be worth to 
focus future Casimir experiments to set bounds on extra-gravitational forces on the actual parallel-plane geometry, 
without the drawbacks of a virtual mapping from the sphere-plane geometry made explicit in this paper. The 
stronger force signal expected for the same distance between the two surfaces, the reduced sensitivity to distance- 
dependent contact potentials due to image charges, the absence of deviations from a uniform radius of curvature, the 
existence of exact mode summation techniques to compute the Casimir force, the possibility to control parallelism 
using recently developed technology |50| , and the possibility to compensate off-line the lack of parallelism by using 
the PFA as discussed in [5l[ , all point in the direction to continue this class of experiments in the actual parallel plane 
configuration, extending below the 10 fim range the results of the experiments described in 0, H, 0, [3, M: [13, [HI, 

Mm 
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